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Abstract 

The physics of organic semiconductors is dominated by the effects of energetic disorder. We 
show that image forces reduce the electrostatic component of the total energetic disorder near 
an interface with a metal electrode. Typically, the variance of energetic disorder is dramatically 
reduced at the first few layers of organic semiconductor molecules adjacent to the metal electrode. 
Implications for charge injection into organic semiconductors are discussed. 

PACS numbers: 73.30.+y, 73.40.Ns 
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The past two decades have been characterized by dramatic advances in the performance 
of organic semiconductor devices, giving rise to the field known as organic electronics Q|. 
Light-emitting diodes P], thin film transistors P], and photovoltaic cells 0|, are examples 
of devices being developed based on organic semiconductors. Critical to the operation of 
all these devices is the process of charge injection from a metal electrode into the organic 
semiconductor. The efficiency of organic light emitting diodes, for example, is directly 
related to the ability of the contacts to supply the organic bulk with charge p . Despite the 
great technological importance of charge injection, the physics of this process remains poorly 
understood. This may be ascribed to the fact that transport in organic semiconductors is 
very different from that in their inorganic counterparts. In the former materials, transport 
takes place by hopping between localized electronic states, which are distributed in energy 
due to spatially correlated energetic disorder with the standard deviation cr ~ 0.1 eV 0,01. 

Recently, it was recognized that energetic disorder in organic materials used in today's 
devices affects the injection efficiency P, 0, Q]. First, disorder was shown to increase 
injection and, second, it was proposed as a major reason for the unusually weak temperature 
dependence of the injected current 11. 1^. 

The injection properties of metal/organic interfaces depend on the properties of a thin 
organic layer directly contacting with the metal. It is well known that the structure of this 
interface layer is typically different from the bulk structure of the organic material. For 
this reason we may suspect that the variance a'j of the energetic disorder at the interface 
differs from the variance al of the disorder in the bulk of the organic material. In literature 
0, ^, Q] , a calculation of the effect of energetic disorder on the injection has been carried 
out using bulk disorder parameters (basically, its variance al). To some extent this could 
be explained by lack of any detailed knowledge of the structure of the interface layer. In 
addition, experimental data of temperature dependence of the injected current seem to 
supports the idea that cr,/ ~ ah Q]. At the same time, it is well known that frequently a 
surface dipole layer is formed directly at the interface, providing an abrupt leap in carrier 
energy in the range of 0.3 — 1 eV It is reasonable to assume that such a layer has some 
degree of disorder and, thus, induces additional energetic disorder in neighboring layers of 
organic materials The magnitude of this energetic disorder should decay while going 
away from the interface, so cXj > 0"^, but this magnitude is completely unknown; in the 
calculations carried out in Ref. very speculative parameters have been used to estimate 
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CTj. 

We are going to show that this problem of the relative magnitude of and at has 
an additional and quite remarkable twist, because in organic devices sandwiched between 
conducting electrodes the bulk disorder itself depends on the proximity to the electrode. A 
well known fact is that a significant part of the total energetic disorder in organic materials 
has an e.ect.ostat.c ong.n. Fo. pel. .nate.a. tMs . dipola.. diso.de. flffQ. wKHe fo. 
non-polar materials it is quadrupolar disorder Our major goal is to demonstrate that 
the electrostatic disorder at the vicinity of metal/organic interface differs from the bulk 
disorder far away from the electrode. 

Indeed, the electrostatic energetic disorder in organic materials is directly proportional to 
the disorder in the spatial distribution of electrostatic potential, generated by randomly situ- 
ated and oriented dipoles or quadrupoles. In organic layers sandwiched between conducting 
electrodes this spatial distribution must obey a boundary condition at the electrode surface: 
at this surface the potential should be a constant. Thus, at the electrode surface there is no 
energetic disorder at all, irrespectively to how disordered is the material in the organic layer. 
This means that the magnitude of the dipolar or quadrupolar disorder increases while going 
away from the interface, asymptotically reaching its bulk value. Here we assume that there 
is no significant increase of the dipolar or quadrupolar disorder directly at the interface (i.e., 
a disordered surface dipole layer is absent). Now we are going to support this general idea 
with the calculation of the variance of the dipolar disorder in the vicinity of a conducting 
electrode. 

Let us consider the simplest model of a rigid disordered polar organic material where 
the randomly oriented (and orientationally uncorrelated) point dipoles occupy the sites of 
a simple cubic lattice with the lattice scale a 0, 14|. We consider the vicinity of a metal 
electrode located at 2; = 0, so the lattice occupies the half-space z > with the first lattice 
layer having distance a/2 from the electrode plane. Charge carrier energy at any site m is 
the sum 

U (rm) = e ^ (p{rra, Tn) (1) 

where 0(r, r„) is the electrostatic potential, generated by the dipole, located at the site n. 
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The variance of the disorder is 

where the angular brackets denote the average over positions and orientations of dipoles, 
and the variable C,n equals to 1 if a dipole is located at the site n and otherwise (note that 
(Um) = 0). A spatial average gives 

iUi) = c5ni + (1 - , (3) 

where c is the fraction of sites occupied by dipoles, and taking into account that the average 
over dipole orientations gives (</> (fVn,'r„)) = 0, we obtain 



J2{<PHr^,rr.)). (4) 



From this point the angular brackets denote the orientational average only. In the case of 
an infinite medium without any electrodes 

nr,rn) = — - — (5) 
e\r - r„| 

where e is the dielectric constant of the medium and pn is the dipole moment. In the case of 
semi-infinite medium bounded by an electrode, a boundary condition = at z = has to 
be imposed (we choose the arbitrary constant to be zero). As a result, the source function 
0(r, r„) includes a contribution from the image dipole = {—Pnx,—Pny,Pnz) located at 
= {xn,yn, —Zn)- an average over dipole orientations gives {pniPnj) = \p'^^ij (with the 
obvious modification for (^PniPnj^), and finally 

2^2^ „ r 1 1 

(6) 



Zn>0 



\rn -A Wn + Z\ 



z -p z 



here the vector z = (0, 0, z) and a does not depend on x and y. The lattice site with r„ = z 
is excluded from the sum ©. Note that cr(0) = 0, as it should be, because if the electrostatic 
potential is a constant for 2 = 0, then there is no electrostatic disorder at this plane, no 
matter how many dipoles occupy the half-space 2 > 0. Far away from the electrode 



2/ \ 2 e p c 
a [oo] = a — 



We can perform an approximate analytical summation in Eq. (jH)) if we provide the 
alternative expression for the source function (p (r, ro), in close analogy to the method, used 
in Ref. The source function for the point dipole, located at tq, is the solution of the 

Poisson equation 

A0 = -— p. V,-„5(f-fo) (8) 

and, hence, is proportional to the gradient of the Green function G (r, ro) of the Laplace 
operator with the zero boundary condition at z = 

An 

(j)(f, fo) = — —p- V^oG (f, fo) . (9) 
Replacing summation with integration in Eq. (j3)) we have 



a^iz) ^ : r ^ / dr' 

z'>0 
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where the Green function has the form 
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VpG(2,r') (10) 



G{r, r') = ^J dke''^''~P'^Gk{z, z'), (11) 

, sinhkz- 
Gk{z, z ) = exp{-kz+), 

z^ = max(2;, z'), z_ = min(2;, z'), 

and p = (x, y) and k are two-dimensional vectors. Performing integration over p ' in Eq. 
flTUI) we obtain 
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Here a cut-off at ~ l/ao with ao — a has been introduced to remove the divergence at 
k oo. This cut-off is equivalent to the exclusion of the site with r„ = z in Eq. ©. We 
did not introduce a similar short range cut-off in the integral over z' in Eq. (jl2p because 
this integral converges and the possible correction does not change the result in a qualitative 
way. To obtain an agreement between the bulk ab in Eq. (fT^ and the corresponding exact 
value for the lattice model in Eq. ((7j) we have to set ao ~ 0.76a [l^]- This choice of ao leads 
to a remarkably good agreement between the approximate analytic expression (fT^ and the 




FIG. 1: Dependence of cr'^{z) on z: the solid line present the result of Eq. (|12() for cq = 0.76a, 
while the squares correspond to the direct summation using Eq. 

result of the direct summation according to Eq. in the whole range of distance from the 
interface (see Fig. [T]). 

Note that for the transport sites situated within a distance of 5 — 6 lattice sites to the 
interface, the amplitude of energetic disorder is significantly decreased in comparison to its 
bulk value. Yet this very thin layer is of crucial importance for injection in organic devices. 
We anticipate that the reduction of a should lead to a stronger temperature dependence for 
the injected current density in comparison with the treatment with a{z) = const, provided 
in the paper by Arkhipov et al. According to this model, the injected current density is 



Joe / dzoexp(-27Zo)u;esc(;2o) / dU'Bo\{U')g[Uo{zQ) - U'] , Uo{zo) = A-- 6^2:0 

(13) 

where Bo\{U) = exp{—U/kT) for U > and Bol(f/) = 1 otherwise, and u'esc(^) is the 
probability for a carrier to avoid the surface recombination 
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Here E is the applied electric field, g{U) is the density of states in the organic material (a 
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Gaussian density of states is usually assumed), 7 is the inverse localization radius, A is the 
interface barrier, and d is the minimal distance, separating the electrode surface and the 
first layer of the organic material. A natural generalization of the Eq. (jl3|) to our case is 
straightforward: we have to let the density of states g depend on z through the Eq. (jl2j) . We 
performed the calculation using parameters provided in Ref. n| for the injection of holes 
from the Ag electrode into poly-dialkoxy-p-phenylene vinylene: A = 0.95 eV, 7 = 0.33 A~^, 
= 5 X 10^ V/cm, ah = 0.11 eV, and d = 12 A (it was assumed that d = a). Note that all 
these parameters were used in Ref. for the comparison between the experimental data 
and Eq. |T!?|l not as fitting parameters, but have been measured independently. The result 
of the calculation is shown in Fig. |21 A transformation of the curve occurs as anticipated: 
because of the smaller a at the interface, the temperature dependence of the injected current 
becomes stronger and does not agree with experimental data anymore. In fact, the agreement 
between the experimental points and the curve, calculated using Eq. (fT^ for a{z) = const, 
is not as perfect as it appears to be, because we have to expect that d < a is a better choice 
for the minimal distance to the electrode (a is the intersite separation). For d < a the curve, 
calculated by Eq. (fT!^ for a{z) = const, goes up (see Fig. |21 the upper solid curve), and the 
agreement worsens. Additionally, small distances to the electrode are very important in the 
integral (fT^ . so the use of the macroscopic e in Eq. (fT^ is dubious. Again, the decrease in 
e moves the curve for a(z) = const in Fig. |21 upwards. 

Q 

All these reservations nothwithstanding, if we believe that the model of Arkhipov et al. 
is valid, then the significant discrepancy between the lower broken line and the experimental 
points in Fig. |21 seems to be an indication of the need to introduce an additional disorder 
(with a ^ 0.1 eV) at the interface. The possible source of this additional disorder could 
be the surface dipole layer. From this point of view, the experimental results, provided in 
Ref. Ill and connected to our consideration, in fact support the idea of a disordered surface 
dipole layer: we clearly need additional disorder at the interface to compensate the decrease 
in the electrostatic disorder, provided by the bulk molecules. 

Possible generalizations (taking into account the possible additional spatial disorder at the 
interface, roughness of the metal/organic interface) do not change our main conclusion that 
the electrostatic part of the energetic disorder in organic materials, provided by molecules 
in the bulk of the material, is significantly suppressed at the interface. 

If a disordered surface dipolar layer is indeed formed at the interface, the picture suggested 
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FIG. 2: Plot of the temperature dependence of the injected current density Jt (normalized by 
the corresponding density J300 for 300 K) for the injection of holes from the Ag electrode into 
poly-dialkoxy-p-phenylene vinylene, calculated using Eq. ((T^ for a{z) = const (solid lines) and 
(t{z), calculate by Eq. (|12() (broken line), correspondingly. The squares indicate the experimental 



points from Ref. 



The upper solid line has been calculated for d = a/2, and the lower one for 



d = a, correspondingly. 

in this paper has to be modified. In this case the magnitude of the total energetic disorder 
could decrease with the increase of the distance to the electrode, or it could still increase, 
depending on the relative amplitudes of the bulk and surface contributions. 

In conclusion, we have shown that the energetic disorder in organic semiconductors de- 
creases dramatically in the neighborhood of a metal electrode. The ramification of this study 
is that disorder parameters derived from bulk measurements do not describe first few layers 
near the metal. As a result, simple models that predict enhancement of charge injection in 
organic semiconductors due to presence of disorder need to be reexamined. 
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03-33067. The research described in this publication was made possible in part by Award 
No. RE2-2524-MO-03 of the U.S. Civilian Research & Development Foundation for the 
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